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Combined  nonlinear  effects  in  two-photon  absorption  chromophores 

at  high  intensities 

R.  L.  Sutherland*3,  D.  G.  McLean3,  M.  C.  Brant3,  J.  E.  Rogers'3,  P.  A.  Fleitzc,  and  A.  M.  Urbasc 
“Science  Applications  International  Corporation,  Dayton,  OH,  USA  45431 
bUES,  Inc.,  Dayton,  OH,  USA  45432 

cAir  Force  Research  Laboratory  (AFRL/MLPJ),  Wright-Patterson  Air  Force  Base,  OH,  USA  45433 

ABSTRACT 

Large  two-photon  and  excited  state  absorption  have  been  reported  in  donor-acceptor-substituted  n-conjugated  mole¬ 
cules.  We  have  performed  detailed  nonlinear  absorption  and  photophysical  measurements  on  a  system  of  AFX  chromo¬ 
phores  and  calculate  the  nonlinear  transmission  based  on  an  effective  three-level  model.  A  numerical  model  that  in¬ 
cludes  far  wing  linear  absorption  has  been  developed  and  compared  with  an  analytical  three-photon  absorption  model. 
The  models  are  in  accordance  and  yield  excellent  agreement  with  experimental  nonlinear  transmission  data  for  0.02-M 
AFX  solutions  up  to  laser  intensities  ~  1-5  GW/cm2.  We  have  extended  our  modeling  efforts  to  include  some  new  ef¬ 
fects  that  may  be  anticipated  in  this  regime,  such  as  stimulated  scattering,  molecular  interactions,  and  saturation.  Effects 
of  chirped  pulses  and  linewidth  of  the  pump  laser  on  stimulated  scattering  are  included.  Self-focusing  and  de-focusing 
are  also  considered.  We  report  on  our  experimental  observations  of  various  materials  and  discuss  results  with  respect  to 
our  extended  theoretical  models. 

Keywords:  two-photon  absorption,  excited  state  absorption,  nonlinear  transmission,  nonlinear  scattering 

1.  INTRODUCTION 

Symmetric  and  asymmetric  electron-donor/acceptor-substituted,  n-conjugated  systems  are  a  major  class  of  enhanced 
two-photon  absorption  (TPA)  materials.12  The  molecular  TPA  cross  section,  a 2,  is  often  characterized  by  nonlinear 
transmittance  (NLT)  experiments,  both  in  the  nanosecond  and  femtosecond  regimes.  Flowever,  the  nanosecond  meas¬ 
urements  in  these  materials  typically  yield  values  of  n2  larger  by  more  than  two  orders  of  magnitude.3,4  Excited  state 
absorption  (ESA)  has  been  postulated  to  play  a  role  in  the  nanosecond  measurements,  and  for  this  reason  the  nonlinear 
parameters  have  been  called  effective  TPA  cross  sections.5,6  Effective  TPA  cross  sections  are  useful  for  screening 
nonlinear  absorption  (NLA)  effectiveness  of  materials  in  the  nanosecond  regime.  Flowever,  understanding  the  physics 
of  nanosecond  NLA  is  necessary  for  designing  improved  materials. 

TPA  followed  by  ESA  was  observed  as  early  as  1974,  and  a  value  for  the  product  of  the  ESA  cross  section  and  excited 
state  lifetime  of  the  chromophore  was  estimated  based  on  a  rate  equation  analysis.7  Several  authors  have  posited  the 
effects  of  excited  states  on  nonlinear  absorption  and  refraction,  and  these  effects  have  been  convincingly  demonstrated 
in  degenerate  four-wave  mixing  and  Z-scan  experiments.811  Evidence  of  ESA  has  been  demonstrated  in  NLT  and  Z- 
scan  measurements  of  organic  materials,  and  effective  (intensity  dependent)  TPA  coefficients  have  been  employed  to 
estimate  ESA  cross  sections.3 11  Recently,  we  have  independently  identified  and  characterized  these  excited  states,  then 
used  this  information  to  theoretically  model  nanosecond  NLT  measurements.12  L'  The  nature  of  these  excited  states  is 
important.  For  example,  properties  of  the  excited  singlet  state  are  significant  for  fluorescence  imaging  and  laser  applica¬ 
tions,  whereas  the  triplet  state  plays  an  important  role  in  photopolymerization  and  photodynamic  therapy.  Moreover,  the 
intrinsic  (pulse-width  independent)  value  of  the  TPA  cross  section  a2  as  well  as  the  significance  of  ESA  in  comparison 
to  other  potential  competing  or  contributing  nonlinear  mechanisms,  such  as  stimulated  scattering,  self- 
focusing/defocusing,  and  quenching,  need  to  be  elucidated  for  accurate  modeling  of  nonlinear  absorption  in  the  nano¬ 
second  regime.  Accounting  for  these  other  phenomena  as  potential  sources  of  NLT  would  enhance  our  understanding  of 
the  process  in  these  chromophores. 
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We  have  performed  detailed  nonlinear  absorption  and  photophysical  measurements  on  a  system  of  AFX  and  similar 
chromophores  and  calculate  the  nonlinear  transmission  based  on  an  effective  three-level  model.  A  numerical  model  that 
includes  far  wing  linear  absorption  has  also  been  developed  and  compared  with  the  analytical  three-photon  absorption 
model.  The  presence  of  small,  linear  absoiption  may  lead  to  a  two-step,  two-photon  absorption  process.14  The  models, 
however,  are  in  accordance  and  yield  close  agreement  with  experimental  nonlinear  transmission  data  for  0.02-M  AFX 
solutions  up  to  laser  intensities  ~  1-5  GW/cm2.  We  have  extended  our  modeling  efforts  to  include  some  new  effects  that 
may  be  anticipated  in  this  and  higher  intensity  regimes,  such  as  stimulated  scattering,  molecular  interactions,  and  self- 
(de)focusing.  We  report  here  on  these  extensions  to  our  models  and  discuss  the  effects  these  additional  phenomena  can 
be  expected  to  have  on  NLT  measurements. 

2.  MATERIALS 

The  chromophores  used  in  this  study  are  part  of  a  class  of  TPA  materials  (designated  AFX)  having  a  design  based  on  a 
push-pull  charge-transfer  model  and  multidimensional  conjugation  motif.6  Examples  of  these  molecules  are  given  in  an 
accompanying  paper.15  These  molecules  were  synthesized  and  purified  in  the  Polymer  Branch  of  the  Air  Force  Materi¬ 
als  Directorate  (AFRL/MLBP)  as  described  in  References  4,  6,  and  22.  Solutions  of  each  were  prepared  in  tetrahydrofu- 
ran  (TF1F)  for  both  photophysical  characterization  and  NLT  experiments. 
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Figure  1.  Linear  absorption  spectrum  (ground  state  cross  section  ctso)  of  AF350  in  THF.  The  circles  are  data  and  the  curve  is  a  Voigt 
function  fit  to  the  data.  Spectra  for  other  AFX  molecules  are  similar. 

A  typical  linear  absorption  spectrum  is  shown  in  Figure  1.  The  /,max  is  in  the  UV  or  deep  blue  region  of  the  visible  spec¬ 
trum,  with  absorption  band  edges  ranging  from  ~  450  nm  to  ~  470  nm.  In  each  case  there  is  very  little  absorption  (cjSo  < 
10"21  cm2)  in  the  2. 5-1. 4  eV  (500-900  nm)  region.  The  purity  should  be  at  least  98%  for  glassy  materials  (e.g.,  AF455) 
and  99%  for  crystalline  compounds  (e.g.,  AF240)  based  on  elemental  analysis  together  with  spectroscopic  methods  and 
melting  point  determination.  There  is  only  a  single  peak  detected  by  liquid  chromatography  for  materials  such  as  AF350 
and  AF455.  An  impurity  with  a  peak  in  the  500-900  nm  range  with  a  peak  absorption  cross  section  of  >  10"17  cm"2 
would  be  observed  at  a  fractional  concentration  of  10  ppm  or  less.  We  conclude  that  it  is  unlikely  that  there  is  signifi¬ 
cant  impurity  absorption  in  the  samples. 
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3.  THEORETICAL  MODELS 


3.1.  Effective  TPA  and  3PA  Models 


Simple  analytical  models  that  can  adequately  explain  and  reliably  predict  the  nanosecond  NLT  based  on  measurable 
properties  of  the  chromophores  and  the  laser  pulse  are  desirable.  The  model  we  have  developed  for  this  purpose  is  simi¬ 
lar  to  that  used  for  reverse  saturable  absorption  media,  but  with  TPA  from  the  ground  state.12  We  consider  ESA  from 
the  lowest  lying  singlet  and  triplet  states.  Two-photon-induced  ESA  is  a  three-photon  process,  although  three  photons 
are  not  absorbed  simultaneously  due  to  the  finite  lifetimes  of  the  excited  states.  In  the  following  we  present  an  analytical 
model  that  incorporates  these  ideas  as  an  effective  three-photon  absorption  (3PA).12 

We  have  for  the  spatial  evolution  of  the  intensity  /, 


^=-(P  +  Yeff/K2 


where 


Tp°eff02^ 


Yeff  2hm 
°eff  =  TFsi  +(j-nW 


(1) 

(2) 

(3) 


P  =  a2N is  the  TPA  coefficient,  and  q  =  (ts/tp){1  -  (ts/tp)[1  -  exp(-ip/TS)]}.  We  have  thus  modeled  the  system  as  a  two- 
photon  resonant,  three-photon  absorption  medium,  with  an  effective  intermediate  state  that  represents  a  time-averaged 
combination  of  the  Si  and  T!  states.  Several  authors  have  taken  the  quantity  in  parentheses  in  Eq.  (1)  to  represent  an 
effective  TPA  coefficient  poff,  from  which  an  effective  TPA  cross  section  can  be  derived.  We  will  consider  both  effec¬ 
tive  TPA  and  3PA  models  in  analyzing  the  experimental  data.  The  quantity  q  is  the  pulse-averaged,  relative  weight  of 
the  Si  state  contribution  to  the  effective  ESA  cross  section.  When  xp  «  is,  q  ~  Vi  and  ESA  is  essentially  from  the  Si 
state  only.  On  the  other  hand,  when  xp  »  is,  q  ~  0  and  all  ESA  is  effectively  from  the  Ti  state.  When  xp  «  Xs,  the 
growth  of  the  Si  population  is  linear,  and  the  average  growth  time  over  a  rectangular  pulse  is  just  xp/2.  The  maximum 
relative  weight  of  'A  for  either  state  reflects  the  fact  that  in  either  time  regime  the  growth  rate  of  the  respective  excited 
state  population  (singlet  or  triplet)  is  approximately  constant.  The  average  population  density  over  the  entire  pulse  is 
then  just  one-half  of  the  maximum  density  obtained  at  the  end  of  the  pulse,  Nsl  max  =  ^r'^Ti.max  =  TP°2^peak/2^®  ■  In 
an  intermediate  regime  where  xT  »  xp  and  xp  ~  Xs,  q  <  'A  means  that  the  Si  population  at  the  end  of  the  pulse  is  smaller 
than  Ns lmax  because  of  continuous  transitions  to  the  S0  and  Ti  states  as  Si  is  being  pumped,  analogous  to  a  leaky  ca¬ 
pacitor. 

Equation  (1)  can  be  solved  analytically,  but  the  resulting  expression  is  transcendental.  However,  if  we  take  P  «  yeg I, 
the  result  reduces  to  a  simple  closed  form  identical  to  that  for  three-photon  absorption.  Assuming  now  Gaussian  spatial 
and  temporal  profiles  for  the  incident  intensity,  we  integrate  the  resulting  expression  over  space  and  time  to  obtain  the 
energy  transmittance  for  a  sample  of  thickness  d\ 


T  = 


'fxpo 


+oo  _ _ 

J  In  yjl  +  pi  exp(-  2x2 )  +  p0  exp(-  x2 ) 

J  -oo  L 


dx 


(4) 


where  T0  is  the  net  linear  transmittance  from  air  into  the  medium  due  to  Fresnel  losses  only  (including  all  dielectric  in¬ 
terfaces),  p0  =  ( 2yeff7o2/o2fi01 12 ,  and  70  is  the  incident  peak,  on-axis  intensity.  The  parameters  determining  yeff  can  be  found 
from  femtosecond,  picosecond,  and  nanosecond  photophysical  measurements  and  standard  TPA  measurements  in  the 
femtosecond  regime. 

For  screening  purposes  it  is  often  useful  to  define  an  effective  TPA  coefficient  as  above.  The  effective  TPA  cross  sec¬ 
tion  can  be  found  by  fitting  the  NLT  data  to  the  following  equation: 
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(5) 


T  = 


where  q0  =  Pc(T7Vo'7.  The  effective  TPA  cross  section  is  then  defined  by  o2eff  =  P IN. 


3.2.  Numerical  ESA  Model 


The  numerical  model  calculates  the  spatial,  time,  and  radial  dependence  of  three  state  populations;  all  of  the  transitions 
between  these  states;  the  bimolecular  processes  of  triplet-triplet  annihilation,  ground  state  quenching,  and  oxygen 
quenching.1'  The  radiation  transport  equation  is  then, 


dl_ 

dz 


—  ^SO^SO^  <TT\NT\I  ■ 


The  population  rate  equations  are 


SNs  i  _  / 
dt  Pico 


rr?n  4” 


&2I 


N 

I Nso  t  KTT  N T [ 


dNjn  _  Vt  jy  _ 
dt  re  Sl 


+  Ktt  Nti  +/<^ST^S0 


I NT 


SNo 

dt 


kot^o^t\  +k0(N0  N0) 


(6) 


(7a) 

(7b) 

(7c) 


N 


so 


N-Nsl-Nn. 


(7d) 


where  ktt,  kSt,  Kot,  N0,  k0  are  the  triplet-triplet  annihilation  rate,  the  ground  state  self  quenching  rate  of  the  triplet  state, 
the  oxygen  quenching  rate  of  the  triplet  state,  the  oxygen  concentration,  and  the  oxygen  singlet  state  relaxation  rate  re¬ 
spectively.  The  rate  equations  are  solved  at  every  time  step  using  a  stiff  Runge-Kutta  algorithm  supplied  in  MathCad. 
These  populations  are  used  to  calculate  the  change  in  the  intensity,  which  is  then  rescaled  assuming  Gaussian  beam 
propagation  at  each  z  step.  A  range  of  input  intensities  are  calculated  and  used  to  calculate  both  the  NLT  variation  and 
the  radial  dependence.  A  Gaussian  time  profile  is  used.  The  radial  and  temporal  integrations  are  performed  to  yield 
transmitted  pulse  energy.  A  normal  reflection  Fresnel  loss  is  imposed  on  the  incident  intensity  and  again  on  the  output 
intensity. 


3.3.  Stimulated  Scattering  Model 

Recently,  He  et  al.  reported  the  observation  of  stimulated  backscattering  in  a  two-photon  absorption  medium,  where  the 
frequency  of  the  stimulated  wave  was  identical  to  the  incident  laser  frequency,  and  the  small-signal  gain  was  quadratic 
in  the  incident  laser  intensity.16 17  They  considered  a  stimulated  thermal  Rayleigh  scattering  model  based  on  TPA- 
enhanced  temperature  and  density  fluctuations,  but  ruled  this  out  due  to  (a)  the  broad  linewidth  of  their  pump  laser, 
which  would  severely  reduce  the  gain  of  the  stimulated  wave,  and  (b)  the  fact  that  the  peak  gain  in  this  theory  is  for  an 
anti-Stokes-shifted  wave,  contrary  to  their  experimental  results.  They  concluded  that  the  stimulated  wave  was  due  to  a 
Bragg  grating  formed  by  the  superposition  of  the  incident  laser  beam  and  an  elastically  (Rayleigh)  backscattered  wave, 
creating  an  index  grating  via  a  TPA-resonance  enhanced  nonlinear  index  coefficient  n2 .  The  stimulated  backscattered 
wave  is  then  due  to  reflection  of  the  laser  from  this  grating. 

Two-beam  coupling  (TBC)  has  a  rich  history  in  nonlinear  optics.  Early  studies  by  Silberberg  and  Bar  Joseph  established 
two  conditions  for  one-way  energy  transfer  between  waves  in  a  Kerr  medium:  1)  the  frequencies  must  be  nondegener¬ 
ate,  and  2)  the  nonlinearity  must  have  a  finite  response  time.18  For  a  positive  Kerr  nonlinearity  energy  transfer  under 
these  conditions  is  always  from  the  high  frequency  to  low  frequency  beam.  Yeh  later  gave  an  exact  solution  for  energy 
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transfer  between  two  co-propagating  beams  including  the  effect  of  linear  absorption.20  In  each  of  these  studies,  both 
beams  had  non-zero  input  values. 

TBC  also  describes  several  stimulated  scattering  phenomena,  such  as  stimulated  Raman  (SRS),  Brillouin  (SBS),  and 
Rayleigh-wing  scattering  (SRWS),  where  one  beam  is  derived  internally  from  the  incident  laser.  For  example,  the 
Stokes  wave  in  these  phenomena  arises  internally  from  scattered  light  and  experiences  gain  at  the  expense  of  the  inci¬ 
dent  laser  beam.  In  each  of  these  stimulated  scattering  effects  the  small-signal  gain  of  the  Stokes  wave  is  proportional  to 
the  laser  intensity,  not  the  square  of  the  intensity.21 

The  model  proposed  by  He  et  al.  would  thus  appear  to  violate  the  two  conditions  previously  established  for  one-way 
energy  transfer  by  TBC  in  a  Kerr  medium:  the  frequencies  of  the  two  beams  are  not  different,  and  the  third  order  sus¬ 
ceptibility  related  to  TP  A,  due  to  an  electronic  polarization,  has  a  virtually  instantaneous  response  for  nanosecond  laser 
pulses. 

In  Ref.  23  Sutherland  considered  a  field  with  a  time  dependent  amplitude  and  phase,  E(t)~  rt(/)exp[- (>(/)].  A  Taylor 
series  expansion  of  the  phase  yields  (ignoring  a  constant  term)  <p(t)=  cot  +  bt2+...,  where  co  is  the  central  frequency  of 
the  wave,  and  b  is  a  linear  chirp  coefficient.  For  simplicity,  higher  order  terms  were  ignored. 

Let  the  total  field  in  a  medium  of  length  d  be  described  by 

E(z,t)=  A L(z,t  +  T)exp\  i[kz -  6Jt - b(t  +  t)2  ]  }  ^ 

+  As  (z,  t  -  r)  exp  j  i[-kz  -  cot -b(t-r )2  ]  ]+  c.c. 


where  k  =  nw/c,  n  is  the  linear  refractive  index,  and  2r=  2 n(d  -  z)!c  is  the  relative  time  delay  at  position  z  and  time  t 
between  the  forward  propagating  laser  wave  (L)  and  the  backward  propagating  scattered  wave  (S).  We  assume  that  the 
scattered  wave  originates  at  z  =  d  by  some  elastic  scattering  process,  so  As(d,t)=  yfrjAL(d,t)  where  77  is  a  constant  «  1. 
Consequently,  the  scattered  wave  has  the  same  spectral  composition  as  the  incident  laser  wave.  However,  for  z  ^  d  a 
lower  frequency  part  of  the  scattered  wave  is  always  interacting  with  a  higher  frequency  part  of  the  incident  wave.  The 
total  polarization  of  the  medium  is  given  by 

P  =  e*\?f  +(NJNW')  +3  X{i){E2)]E  (9) 

where  /^")  is  the  n-th  order  susceptibility,  and  the  angular  brackets  indicate  an  average  over  a  time  longer  than  an  optical 
period  but  shorter  than  (2 bz)A .  We  have  assumed  that  TPA  produces  a  single  excited  state  of  number  density  Ne  «  N, 
the  total  number  density  of  the  nonlinear  chromophore.  Both  A%fl)  =  ~Zg>  and  are  complex  quantities  ( g  signi¬ 
fies  the  ground  state  of  the  medium).  The  population  of  the  state  excited  by  TPA  decays  back  to  the  ground  state  with  a 
time  constant  Te,  so  Ne  obeys  the  following  kinetic  equation: 


dNe  _a2NI2  Ne 
dt  2Tm  Te 


(10) 


where  l(z,t)  =  2eanc[E2(z,tfj  is  the  total  intensity.  Let  the  amplitudes  be  slowly  varying  in  time  compared  to  Te.  Equa¬ 
tion  (10)  can  then  be  integrated  to  yield 

Ne/N  =  C0  +{c,  exp[i(2fe-4Zm)]  +  C2  exp[i(4fo-8Z>rt)]+c.c.} ,  (Ha) 


C0  = 


°lTe 

2ha> 


(^i  +^s  +^bLIs), 


C,  = 


<?2Te 

2ha> 


4gp ncALA*s  N 
\—iAbxTe 


a2Te  (2 s0ncALAsf 
2  hco  \-i%brTe 


(lib) 

(11c) 

(lid) 
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where  So  is  the  free-space  permittivity.  We  will  make  the  assumption  that  ( AbrTe )2  «  1. 

Following  a  procedure  directly  analogous  to  that  of  Yeh20  and  Boyd21  for  non-degenerate  TBC,  Eqs.  (lla)-(lld)  are 
substituted  into  Eq.  (9),  and  then  both  Eqs.  (8)  and  (9)  are  substituted  into  Maxwell’s  wave  equation  in  the  slowly  vary¬ 
ing  amplitude  approximation.  Matching  up  synchronous  terms,  the  following  coupled-wave  equations  for  the  laser  and 
backscattered  intensities  are  derived: 


=  -g(l  -  z/d\lL  +  Is  )ILIS  -  yeff{l2L  +  31  $  +  6 ILIS  ]lL  -  p{lL  +  21  s  )IL 
az 

~~r~  =  _g(l —  z/d  \li  +  Is  )IlIs  +  Yeff  (-^z  +Is  +  ^dLIs  \s  +  fifai.  + 1 s  Vs 


(12a) 

(12b) 


where  g  and  arc  the  backward  wave  gain  and  effective  three-photon  absorption  (3PA)  coefficients,  respectively,  with 


g  = 


8bTea>A%<jPd 


(13) 


NAo 

I eff  ~  ~~p.  ’ 


(14) 


where  Axf  =  Rc(a^(1)),  Act  =  ae-<r„ 
ground  states,  and  Imt  =  (2ha>/cr^Te)'/~ 
(2). 


is  the  difference  between  the  linear  absorption  cross  sections  of  the  excited  and 
is  the  two-photon  saturation  intensity.  Note  the  similarity  between  Eqs.  (14)  and 


Linear  chirp  is  not  a  unique  requirement  for  the  energy  transfer  between  incident  laser  and  elastically  backscattered 
waves.  Another  potential  mechanism  involves  multimode  beams.  Take  the  simple  case  where  the  incident  laser  wave 
consists  of  two  modes:  coi,  and  a>i  =  co0  +  Sco>  a>0  ( So)  «  co 0).  Flere  the  modes  are  not  in  any  time -ordered  sequence 
(such  as  a  chirp,  where  the  frequency  varies  linearly  in  time),  and  can  thus  interact  simultaneously  through  the  nonlinear 
polarization  of  the  medium  over  the  duration  of  the  pulse.  Employing  the  mechanism  involving  a  TPA-populated  ex¬ 
cited  state  described  above,  there  will  be  energy  transfer  from  the  a>\  mode  of  the  incident  laser  beam  to  the  (ot)  mode  of 
the  backscattered  beam.21  Likewise,  the  a>\  mode  of  the  backscattered  wave  will  yield  its  energy  to  the  (0„  mode  of  the 
laser  wave.  Thus,  we  need  consider  only  coupling  between  the  laser  frequency  at  a>i  with  the  backscattered  wave  at  a>0. 
In  this  case  the  above  formalism  may  be  used  with  the  substitution  4b  r— >  Sco. 


3.4.  Nonlinear  Beam  Propagation  Model 

To  incorporate  the  effects  of  nonlinear  focusing  on  the  NLT,  we  must  solve  the  foil  three-dimensional  wave  equation  in 
the  slowly  varying  amplitude  approximation.  The  solution  for  the  amplitude  after  propagation  through  a  small  distance 
Az  along  the  z-axis  can  be  written  formally  as24,25 

,  z  +  Az)  =  expj^'  H^A{r,z)  (15) 

where  H  is  an  operator  containing  radial  and  nonlinear  dependence  of  the  beam  propagation.  We  assume  the  beam  has 
axial  symmetry.  Thus, 


^  +  (16) 
or  ~  r  or 

The  function  /  contains  the  nonlinear  terms  of  the  wave  equation.  In  an  isotropic  system,  this  will  include  only  odd- 
order  nonlinear  susceptibility  terms.  To  third  order,  it  will  include  the  effects  of  nonlinear  refraction  and  TPA.  We  wish 
to  also  incorporate  the  effect  of  excited  state  absorption  through  the  effective  3PA  approach.  Hence,  the  nonlinear  func¬ 
tion  in  Eq.  (16)  becomes 
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(17) 


/  = 


%nn2n2  .  2mB  4 ny  „ 
- —^  + 1 - — +  / - —  AA 


A 


nA 


AA 


where  «2  is  the  nonlinear  refractive  index  coefficient  in  intensity  units  (cm2/W). 

We  used  the  technique  of  Feit  and  Fleck24  to  solve  the  wave  equation,  which  involves  expanding  the  exponential  opera¬ 
tor  of  Eq.  (15)  in  a  Taylor  series  and  applying  the  operators  term  by  term  on  the  wave  amplitude  A(r,z).  The  radial  de¬ 
rivatives  are  found  by  writing  the  wave  amplitude  as  a  digital  Fourier  series  and  using  fast  Fourier  transform  algorithms. 
The  input  amplitude  is  taken  to  be  a  focused  Gaussian  beam  with  the  waist  in  the  center  of  the  sample.  For  free-space 
propagation,  we  kept  terms  to  second  order  in  the  expansion  of  Eq.  (15).  We  checked  the  accuracy  of  this  approach  by 
computing  the  linear  propagation  of  a  focused  Gaussian  beam  and  comparing  it  to  the  analytical  result.  Agreement  was 
better  than  1%  when  including  terms  to  second  order. 


4.  EXPERIMENTAL 

Micromolar  solutions  were  prepared  for  conventional  photophysical  measurements  (lifetimes,  excited  state  cross  sec¬ 
tions,  and  quantum  yields).  Details  of  these  experiments  can  be  found  elsewhere.26  All  solution  samples  for  nanosecond 
NLT  measurements  had  concentrations  of  0.02  M  (mol/L)  and  were  placed  in  1-mm  glass  or  fused  silica  cuvettes.  In¬ 
trinsic  02  values  were  obtained  from  independent  femtosecond  measurements.12 

Nanosecond  nonlinear  transmittance  measurements  were  performed  with  a  Nd:YAG  pumped  optical  parametric  oscilla¬ 
tor  (OPO)  tuned  from  660  to  880  nm.  The  pulse  was  Gaussian  shaped  with  tl  =  3.5  ns.  The  beam  was  focused  with  an  f 
=  50-cm  lens  into  the  sample.  Over  the  length  of  the  sample  (1  mm)  the  beam  was  essentially  collimated.  The  beam 
shape  was  slightly  elliptical,  with  a  geometric-mean  1/e2  radius  w  ~  18.4  -  18.9  pm,  assuming  an  approximately  Gaus¬ 
sian  beam  shape.  The  energy  was  varied,  and  incident  and  transmitted  energies  were  measured  with  energy  meters.  To 
rule  out  the  effects  of  output  beam  spreading  on  transmittance  measurements,  a  large-area  (~  1  cm2)  detector  was  placed 
near  the  exit  of  the  sample  to  collect  all  transmitted  energy.  We  also  looked  for  stimulated  backscattering  by  rotating  the 
sample  slightly  to  avoid  Fresnel  reflected  light  and  measuring  180°-scattered  light  with  an  energy  meter. 

Intermolecular  quenching  properties  were  examined  for  AF350  and  AF455.  The  ground  state  concentration  was  varied 
from  10  pM  to  20  mM,  and  the  triplet  state  lifetime  measured.  At  ~  10  pM  the  lifetime  of  AF350  was  182  ps  under  de- 
oxygenated  conditions,  and  the  lifetime  was  96  ps  at  20  mM.  A  bimolecular  self-quenching  rate  constant  of  2  x  102 
L/mol-s  was  measured.  In  addition,  we  have  measured  the  triplet-triplet  quenching  rate  in  AF455  to  be  ktt  =  2.1  x  109 
L/mol-s.  Thus  it  appears  that  triplet-triplet  annihilation  is  much  more  significant  than  self-quenching  in  these  systems. 
Flowever,  the  triplet  yield  in  these  systems  is  low.  Based  on  an  analysis  of  population  kinetics,  we  would  expect  that  Nj 
<  0.1  A  even  at  the  highest  input  energies.  This  would  yield  a  rate  ktt/Vt  <  4  x  106  s"1  as  an  upper  limit  for  triplet 
quenching,  which  is  of  the  order  of  the  phosphorescent  rate.  Since  this  implies  a  triplet  lifetime  long  compared  to  the 
laser  pulse,  and  since  the  dominant  contribution  to  ESA  appears  to  be  from  the  Si  state,  we  conclude  that  triplet  quench¬ 
ing  processes  have  no  significant  effect  on  the  measured  NLT  in  these  AFX  systems. 


5.  RESULTS  AND  DISCUSSION 


5.1.  TP  A/ESA 

A  common  method  in  the  literature  for  screening  TPA  materials  in  the  nanosecond  regime  is  to  fit  the  NLT  data  to  a 
TPA  transmittance  (see  Eq.  (5))  and  extract  an  effective  TPA  cross  section.  We  give  the  results  of  such  fits  for  two  se¬ 
ries  of  dipolar  (AF240  and  AF270),  quadrupolar  (AF287  and  AF295),  and  octupolar  (AF350  and  AF380)  molecules  in 
Table  1.  Effective  02  as  a  function  of  wavelength  for  a  related  chromophore  (El-BTF)  is  shown  in  Figure  2.  El-BTF  is 
related  to  AF287,  but  has  a  platinum  core  bridging  the  two  benzothiolazole  fluorine  arms.15  Data  such  as  these  are  use¬ 
ful  for  observing  trends  in  molecular  structure,  but  do  not  yield  sufficient  information  to  give  directions  for  further  im¬ 
provement.  For  example,  a  large  effective  02  could  be  due  to  a  large  intrinsic  a2,  a  large  singlet  and/or  triplet  cross  sec¬ 
tion,  and/or  a  large  triplet  yield,  or  perhaps  some  other  not  yet  recognized  nonlinear  effect.  Flence,  these  NLT  measure¬ 
ments  must  be  supplemented  with  photophysical  measurements.  Both  effective  3PA  and  numerical  ESA  models  may 
then  be  applied  to  ascertain  the  quality  of  the  model  in  reproducing  the  NLT  results. 
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Table  1.  Effective  TPA  Cross  Sections  of  a  Series  of  AFX  Chromophores 


Material 

Effective  o2 
(10'2°  cm4/GW) 

AF240 

50 

AF287 

99 

AF380 

114 

AF270 

29 

AF295 

78 

AF350 

139 

Wavelength  (nm) 

Figure  2.  Effective  TPA  cross  sections  of  El-BTF  measured  with  nanosecond  laser  pulses. 

Photophysical  data  at  800  nm  for  this  set  of  AFX  chromophores,  as  well  as  data  for  AF455  at  various  wavelengths,  are 
presented  in  Table  2.  Using  these  data,  and  no  adjustable  parameters,  we  compute  the  expected  transmittance  of  solu¬ 
tions  in  TF1F,  based  on  Eqs.  (2)-(4)  and  compare  these  to  experimental  nanosecond  NLT  data  in  Figure  3.  For  all  of  the 
solutions  in  these  plots,  the  concentration  is  nominally  0.02  M  with  the  exception  of  AF287,  which  is  approximately 
0.01  M  due  to  solubility  limitations.  We  can  immediately  see  that,  although  the  model  predicts  the  general  trends  in  the 
general  pulse  energy  regime,  the  quantitative  agreement  with  the  data  is  not  good  in  all  cases.  This  is  particularly  sur¬ 
prising  for  AF240  and  AF350  for  which  we  obtained  much  better  agreement  between  theory  and  data  at  820  nm.  Before 
addressing  other  possible  nonlinearities  contributing  to  this  discrepancy,  we  should  point  out  that  for  each  of  these 
molecules  the  triplet  yield  is  very  low,  and  thus  the  excited  singlet-singlet  absorption  is  playing  a  major  role  in  the  ESA. 
Excited  state  cross  sections  in  the  800  nm  regime  are  difficult  to  measure  because  this  is  the  fundamental  wavelength  of 
the  Ti:sapphire  laser  used  in  the  experiments.  A  cutoff  filter  is  used  to  eliminate  this  signal,  making  the  measured  ab¬ 
sorbance  signals  at  this  wavelength  somewhat  noisy.  The  signals  for  shorter  wavelengths  do  not  have  this  problem.  For 
example,  we  show  in  Figure  4  transmittance  data  and  model  results  for  AF455  at  wavelengths  of  740,  780,  and  800  nm. 
The  agreement  between  theory  and  data  is  very  good.  We  conclude  that  the  accurate  measurement  of  the  excited  state 
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cross  section  Csi  is  crucial  for  these  AFX  chromophores,  and  measurements  near  the  fundamental  wavelength  of  the 
laser  (800  nm)  need  to  be  improved. 


10'2  10'1  10°  101  102  10'1  10°  101  102  103 


Energy  (pJ) 

Figure  3.  Nanosecond  NLT  data  and  model  results  for  six  AFX  chromophores. 

A  closer  examination  of  the  plots  in  Figures  3  and  4  reveals  that  there  are  also  some  differences  in  the  shape  of  the 
transmittance  curves  between  theory  and  data.  This  is  not  so  surprising  at  the  highest  pulse  energies  where  it  is  conceiv¬ 
able  that  higher  order  nonlinearities  may  be  evident,  but  there  is  also  a  consistent  overestimation  of  the  transmittance  by 
the  theory  in  the  lower  energy  region  where  this  measurement  just  begins  to  roll  over  with  energy.  Although  this  dis¬ 
crepancy  is  reasonably  small  and  somewhat  amplified  due  to  the  semi-log  nature  of  the  plot,  it  appears  to  be  systematic 
and  may  suggest  the  presence  of  other  nonlinearities  that  could  possibly  play  a  significant  role  as  the  pulse  energy  in¬ 
creases.  One  could  question  the  nature  of  the  approximations  leading  to  the  analytical  effective  3PA  model,  but  the 
close  agreement  between  this  analytical  result  and  a  finite  difference  numerical  calculation  of  TPA  followed  by  ESA 
(i.e.,  including  detailed  time -dependent  behavior  of  singlet  and  triplet  populations  and  pulse  intensity  absorption)  indi¬ 
cates  that  the  analytical  approach  is  quite  adequate  for  modeling  these  nonlinear  effects.12  Independent  verification  of 
this  agreement  with  a  generalized  theoretical  treatment  of  laser  pulses  interacting  with  multi-photon  absorbers  strength- 
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ens  this  conclusion.27  We  note  that  a  two-step  TPA  mechanism  (similar  to  reverse  saturable  absorption)  has  been  pro¬ 
posed  to  explain  the  NLT  in  these  types  of  chromophores,14  and  we  have  performed  sensitive  measurements  of  the  lin¬ 
ear  absorption  cross  section  in  AF240,  AF350,  and  AF455.  Including  the  effects  of  these  small  measured  ground  state 
cross  sections  on  the  nonlinear  absorption  (and  population  kinetics),  however,  makes  only  very  small  corrections  to  the 
effective  3PA  theory13  and  does  not  account  for  the  discrepancies  observed  in  the  data.  We  conclude  that  the  analytical 
effective  3PA  result  yields  excellent  quantitative  predictions  for  the  specific  model  in  question,  which  is  TPA  followed 
by  ESA  from  the  singlet  and  triplet  manifolds  of  the  molecules.  Therefore,  we  must  look  at  potential  alternative  proc¬ 
esses  that  may  be  at  work  in  conjunction  with  TP  A/ESA. 


Table  2.  Photophysical  Properties  of  Various  AFX  Chromophores 


Chromophore 

Wavelength  (nm) 

o2  (10-20  cm4/GW) 

CTsi  (10"17  cm2) 

CTti  (10'17  cm2) 

<Pt 

Ts  (ns) 

AF240 

800 

0.19 

1.45 

7.45 

0.064 

2.03 

AF287 

800 

0.50 

4.39 

15.9 

0.074 

1.88 

AF380 

800 

0.63 

6.08 

17.2 

0.130 

1.85 

AF270 

800 

0.15 

4.93 

15.3 

0.031 

1.89 

AF295 

800 

0.32 

5.80 

16.6 

0.050 

1.88 

AF350 

800 

0.49 

6.91 

22.0 

0.035 

1.99 

AF455 

740 

0.37 

1.91 

9.98 

0.060 

2.72 

AF455 

780 

0.81 

1.87 

14.7 

0.060 

2.72 

AF455 

800 

0.51 

1.45 

18.1 

0.060 

2.72 

5.2.  Stimulated  Scattering 

It  is  interesting  that  an  effective  TPA  theory  fits  the  data  better  in  this  low  energy  regime,  indicating  that  some  mecha¬ 
nism  that  has  a  quadratic  intensity  dependence  may  be  at  play.  He  et  al.17  have  observed  this  kind  of  dependence  and 
have  verified  in  their  case  that  this  is  due  to  a  stimulated  scattering  phenomenon.  We  consider  this  next. 

Let  us  examine  Eqs.  (12a)  and  (12b)  for  nonlinear  scattering  in  the  case  when  h<<h  «  constant.  It  can  be  seen  that  the 
backward  stimulated  wave  will  experience  exponential  growth  when  I L  >  /(^g  .  This  defines  the  threshold 

condition  for  stimulated  backscattering  when  there  is  no  linear  absorption.  From  Eq.  (12b),  the  small  signal  gain  G, 
given  in  terms  of  A Is  =  Is(o)-Is(d),  is  G  =  A/s / 1 s(d)r-t-  ij  for  input  intensities  where  the  gain  dominates  TPA.  Also, 
the  reflectance  of  the  scattered  wave  is  defined  by  R  =  Is  (0 )/IL  (o) ,  and  the  change 
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AR  =  R-t7  =  AIs/Il  ccIs(d)lL  oc  GIs(d)/IL  .  These  results  are  in  agreement  with  the  data  presented  by  He  et  al.  for  the 
measured  small-signal  gain  and  reflectance  in  a  0.01-M  solution  of  PRL  802  in  THF.17 


Under  conditions  where  the  nonlinear  absorption  terms  (J3  and  yefj)  can  be  ignored,  i.e.,  where  these  loss  terms  are  small 
compared  to  the  TBC  gain,  Eqs.  (12a)  and  (12b)  can  be  solved  analytically.  The  result  can  be  expressed  as 


R 

(1_/?Xl-/?  +  3/7o)+2/7o2 

1/2 

l 

6l-R  +  2?702 

(l  +  7?)2  exp(r) 

(18) 


where  r]0  =  I^d)/IL( 0),  and 

r  =  Ui-Rfgil(o)d. 


(19) 


Note  that  in  general  r]0  ^  ij,  although  when  the  gain  is  sufficiently  small  tj0  «  i).  Also,  when  the  gain  is  not  too  large  so 
that  (1  -  R)  »  r/0,  Eq.  (18)  can  be  simplified  to  the  following  approximation: 


*M) 

(l  +  T?)2  exp(r) 


(20) 


It  is  interesting  to  compare  this  result  to  the  case  of  SBS,  for  which  the  denominator  of  Eq.  (20)  becomes  exp(r)-// 
with  T  -»(l  -R)gBIL(o)d  ,  where  gB  is  the  Brillouin  gain  coefficient.21  Thus,  in  the  case  of  negligible  loss  by  absorption, 
the  backscattered  reflectance  as  a  function  of  incident  laser  intensity  will  look  similar  to  the  Brillouin  reflectance.  A  plot 
of  the  reflectance  given  by  Eq.  (20)  for  a  constant  input  I^dyhi 0),  over  a  range  of  E  where  the  approximation  is  valid, 
is  shown  in  Figure  5. 


Figure  5.  Backscattered  reflectance  vs.  exponential  gain  factor. 
P  =  0,  yeff=  0  (solid);  g  =  20p2rf,  yeff  =  0  (long  dash);  g=  20P 2d, 
g  =  1 5  Yeff  (short  dash).  Circles  give  R  calculated  by  analytical 
approximation. 


Figure  6.  Nonlinear  transmittance  vs.  incident  intensity. 
Pure  TP  A,  p  =  9.46  cm/GW  (dashed);  p  =  9.46  cm/GW, 
g  =  1800  cm3/ GW2,  Yeff =  120  cm2/GW2,  q  =  0.04  (solid). 
Circles  are  data  according  to  He  et  al.17 


Numerical  solutions  of  Eqs.  (12a)  and  (12b)  are  also  given  in  Figure  5.  These  illustrate  the  agreement  of  the  approxima¬ 
tion  given  by  Eq.  (20)  with  the  exact  result,  and  the  effect  of  including  the  nonlinear  absorption  terms.  Figure  6  gives 
the  laser  transmittance.  Parameters  have  been  adjusted  to  yield  results  close  to  the  experimental  data  of  He  et  al. 16,17  The 
coefficients  used  are  /?  =  9.46  cm/GW  (the  value  quoted  in  Ref.  17),  ij  =  0.04,  g  =  1800  cm3/GW2,  and  yeff=  120 
cm3/GW2,  with  d  =  1  cm.  Here  the  transmittance  is  the  ratio  of  output  power  to  input  power.  To  compute  this,  we  as¬ 
sumed  a  Gaussian  dependence  for  the  input  intensity  and  then  integrated  the  output  intensity  over  the  area  of  the  cylin- 
drically  symmetric  beam.  The  result  is  compared  with  what  would  be  expected  for  pure  TPA  (no  TBC  or  ESA).  The 
departure  from  pure  TPA  in  this  case  is  due  primarily  to  energy  transfer  from  the  laser  to  the  backscattered  beam.  A 
comparison  with  the  data  given  in  Ref.  17  shows  good  agreement  with  the  results  of  this  model. 
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To  get  an  order  of  magnitude  of  the  numbers  involved,  consider  the  experiment  of  He  et  al. 16,17  and  the  results  shown  in 
Figure  6.  For  a  laser  wavelength  of  532  nm  and  Te  ~  1  ns,5,26  Isat  ~  0.7  GW/cm2.  For  a  severely  chirped  pulse  the  laser 
linewidth  is  A®l  ~  2 brL  where  rL  is  the  laser  pulse  width.  For  tsa\J2n~  24  GHz  (0.8  cm"1),  rL  ~  10  ns,  and  N=  6  x  1018 
cm'3  (0.01-M  concentration),  Eqs.  (13)  and  (14)  yield  Aj®  ~  4  x  10"3  and  Acr~  1  x  10"17  cm2.  These  values  yield  ap¬ 
proximate  agreement  with  typical  photophysical  data  for  TPA  chromophores,  and  thus  suggests  that  this  stimulated 
scattering  with  a  quadratic  intensity  dependence  may  be  contributing  to  the  overall  shape  of  the  NLT  curve. 

Notice  that  Rc(^f(3)),  or  n2,  does  not  appear  in  the  gain  coefficient  of  Eq.  (13).  The  reason  for  this  is  that  ^  was  as¬ 
sumed  to  have  an  instantaneous  response  [  [m(j(3))oc  p  ].  There  is  no  two-beam  coupling  in  a  Kerr  medium  with  an  in¬ 
stantaneous  response.20,21  It  is  possible  however,  under  TPA  resonant  conditions,  that  n2  could  have  a  finite  response 
time  (i.e.,  /j/3:i  is  complex  with  a  finite  damping  coefficient21).  However,  if  this  is  included  in  the  development  of  the 
coupled-wave  intensity  equations,  it  would  lead  to  a  term  in  the  gain  G  that  is  proportional  to  IL?  not  /2 ,  which  would 
be  inconsistent  with  the  experimental  results  of  He  et  al.17  It  is  also  quite  likely  that  the  Kerr  refractive  term  would  be 
much  smaller  than  the  term  proportional  to  . 

As  noted  earlier,  the  same  theory  just  described  can  be  applied  to  multimode  pulses  instead  of  chirped  pulses  with  the 
simple  substitution  4br— >  8m.  Pulses  described  in  Ref.  17  appear  to  be  multimode.  However,  as  argued  earlier  the  se¬ 
lective  sign  of  the  TBC  gain  would  then  cause  only  one  mode  of  the  backscattered  wave  to  be  amplified  (assuming  the 
laser  oscillates  on  only  two  modes),  thus  making  the  backscattered  pulse  single  mode.  This  appears  not  to  be  the  case 
experimentally.17  It  seems  likely  that  the  laser  pulse  can  be  both  multimode  and  chirped  and  thus  produce  a  backscat¬ 
tered  wave  in  accordance  with  experiment. 

Applying  this  model  to  our  experimental  data  appears  to  lead  to  a  null  result.  We  have  looked  for  stimulated  backscat- 
tering  by  rotating  the  sample  slightly  to  avoid  Fresnel  reflected  light  and  measuring  180°-scattered  light  with  an  energy 
meter.  For  example,  at  an  incident  pulse  energy  of  124  pj  on  a  sample  of  0.02-M  AF455  in  THF,  we  measured  ~  2  nJ  of 
backscattered  light.  This  would  correspond  to  <  2  x  10"5  efficiency  for  stimulated  scattering,  which  can  be  considered 
negligible  compared  to  the  nonlinear  absorption  loss.  To  date,  we  have  not  observed  any  significant  backscattered  signal 
in  our  NLT  experiments.  The  reasons  for  this  are  probably  related  to  the  nature  of  the  input  pulses.  Whereas  Ref.  17 
used  a  Q-switched,  frequency-doubled  Nd:YAG  laser  as  the  input,  we  use  a  tunable  OPO.  The  OPO  is  pumped  by  a 
seeded,  single  longitudinal  mode  Nd:YAG  laser,  and  the  OPO  pulse  shows  no  evidence  of  multimode  behavior.  Al¬ 
though  the  pump  laser  pulse  may  have  some  chirp  due  to  the  Q-switch,  the  bandwidth  of  the  OPO  pulse  is  dominated  by 
angle  phase  matching  effects.  In  fact,  the  bandwidth  of  the  pulse  is  so  large  that  the  pulse  should  have  a  small  coherence 
length.  Indeed,  measurements  of  the  bandwidth  at  800  nm  (~  37  cm'1)  indicate  a  theoretical  coherence  length  of  only  ~ 
270  pm.  This  is  in  rough  agreement  with  the  loss  of  fringe  visibility  in  our  experiments  when  the  beam  is  split  and 
propagated  over  a  path  difference  of  only  ~  100-200  pm.  This  is  small  compared  to  our  sample  thickness  of  1  mm,  and 
hence  implies  that  the  beam  cannot  interact  via  a  Bragg  grating  coherently  over  a  distance  equal  to  the  sample  thickness. 
We  conclude  that  stimulated  backscattering,  of  the  type  observed  by  He  et  al.  where  the  interaction  length  is  a  factor  ~ 
lOOx  larger  than  our  coherence  length,  is  probably  not  making  a  substantial  contribution  to  the  NLT  measurements  in 
our  experiments  due  to  the  nature  of  the  OPO  pulse. 

5.3.  Nonlinear  Refraction  Effects 

The  usual  assumption  in  the  modeling  of  nonlinear  absorption  is  that  the  beam  does  not  diffract  significantly  during  its 
propagation  through  the  sample.  In  other  words,  we  assume  that  the  only  thing  that  can  change  the  local  intensity  is  ab¬ 
sorption.  However,  since  nonlinear  absorption  is  intensity  dependent,  any  nonlinear  process  that  could  affect  the  local 
intensity  will  in  turn  modify  the  nonlinear  absorption.  In  a  NLT  experiment,  we  normally  design  the  experiment  so  that 
self-(de)focusing  effects  do  not  lead  to  any  loss  in  the  transmitted  (far  field)  beam.  Nonlinear  refraction  can,  however, 
also  affect  the  local  intensity  in  the  sample  and  thus  influence  the  nonlinear  absorption.  To  study  these  effects,  we  have 
employed  the  nonlinear  beam  propagation  model  which  includes  diffraction  and  nonlinear  refraction. 

We  consider  a  hypothetical  material  that  has  all  of  the  photophysical  properties  of  AF455,  but  for  which  we  may  vary 
the  nonlinear  refractive  index  coefficient,  n2.  As  described  in  Section  3.4,  we  include  both  TPA  and  the  effective  3PA 
model.  We  take  the  input  wave  to  the  sample  to  be  a  focused  Gaussian  beam  which  has  its  waist  (for  the  linear  regime) 
located  at  the  center  of  the  sample.  The  Rayleigh  range  of  the  beam  is  longer  than  half  the  sample  length.  For  the  spe- 
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cific  example  considered  here,  we  take  the  beam  waist  to  have  a  radius  w  =  18  pm.  The  wavelength  is  800  nm,  and  the 
sample  is  1  mm  thick.  For  purposes  of  calculating  radial  derivatives  using  the  fast  Fourier  transform  technique,  we  di¬ 
vided  the  beam  radially  into  128  equal  increments  and  selected  a  total  radial  range  of  78  pm.  The  integration  step  along 
the  beam  axis  was  chosen  to  be  Az  =  0.052  =  0.04  pm.  The  output  intensity  was  computed  from  the  output  field  ampli¬ 
tude  and  then  integrated  radially,  assuming  axial  symmetry,  to  compute  the  total  integrated  optical  power. 

Results  of  the  calculation  are  shown  in  Figure  7,  which  plots  the  integrated  power  transmittance  as  a  function  of  inci¬ 
dent  on-axis  intensity  (i.e.,  the  peak  intensity  the  beam  would  have  at  focus  in  a  linear  medium).  The  solid  curve  is  the 
result  of  the  effective  3PA  model,  which  assumes  that  the  beam,  in  the  linear  regime  (low  power),  has  a  constant  radius 
throughout  the  sample.  The  circles  are  the  results  of  the  nonlinear  beam  propagation  model  for  n2  =  0,  i.e.,  negligible 
nonlinear  refraction.  Note  that  at  high  intensities  these  points  he  just  slightly  above  the  curve  for  the  analytical  model. 
This  is  not  an  error  due  to  the  numerical  calculations,  but  is  due  to  the  fact  that  the  beam  propagation  model  includes  the 
effects  of  diffraction:  the  on-axis  intensity  is  slightly  reduced  on  either  side  of  the  focal  region  due  to  beam  spreading, 
whereas  the  analytical  model  assumes  constant  beam  radius.  The  squares  give  the  NLT  that  would  be  measured  if  the 
sample  has  an  n2  equal  to  that  of  carbon  disulfide.  Self-focusing  obviously  increases  the  nonlinear  absorption,  but  the 
variation  is  small  and  probably  within  typical  measurement  uncertainty.  On  the  other  hand,  a  negative  index  such  as 
might  be  associated  with  a  thermal  nonlinearity,  has  a  definite  impact  on  the  transmittance  at  higher  intensity.  Since  the 
medium  becomes  more  absorbing  at  high  intensities,  some  of  this  energy  is  likely  to  be  dissipated  as  heat  and  cause  the 
beam  to  be  thermally  bloomed  to  some  degree.  This  would  increase  the  NLT  and  make  it  appear  that  the  sample  has 
smaller  TPA  and/or  ESA  cross  sections.  A  rather  large  positive  n2  can  make  the  NLT  depart  significantly  from  the  ana- 


Peak  Intensity  (GW/cm2) 

Figure  7.  Nonlinear  power  transmittance  vs.  peak  intensity  for  a  hypothetical  material  like  AF455  but  with  various  values  of  nonlin¬ 
ear  index  coefficient. 

lytical  prediction  even  at  relatively  low  intensity,  as  indicated  by  the  single  diamond  point  in  Figure  7.  It  seems  unlikely 
that  a  third-order  effect  would  contribute  such  a  large  nonlinear  index,  since  n2  is  probably  due  mostly  to  the  solvent  at 
this  concentration,  and  the  n2  for  TFIF  is  smaller  than  that  of  CS2.  Flowever,  ESA  can  make  a  substantial  contribution  to 
the  nonlinear  refractive  index.  From  the  example  given  in  Section  5.2,  the  index  change  due  to  population  of  the  excited 
state  would  be  An  ~  10"\  This  is  a  substantial  index  modulation  and  larger  than  the  one  used  in  the  calculation  in  Figure 
7  (A/7  =  n2I  =  2.6  x  10’4)  for  which  n2  =  8.8  x  10"4  cm2/GW  was  assumed.  Thus,  it  is  conceivable  that  nonlinear  refrac¬ 
tion  due  to  population  redistribution  (i.e.,  the  creation  of  excited  state  species)  may  be  sufficient  to  lower  the  NLT  in  the 
small  energy  regime,  making  it  appear  that  the  medium  has  larger  nonlinear  absorption  cross  sections.  Finally,  we  must 
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also  consider  the  possibility  that  both  positive  and  negative  index  effects  (e.g.,  population  redistribution  and  thermal 
effects)  may  be  at  play  simultaneously  (with  possibly  different  relative  strengths  at  different  pulse  energies),  affecting 
the  overall  shape  of  the  NLT  curve.  This  is  only  speculative  at  this  point,  but  plausible  based  on  the  above  analyses. 
More  work  is  necessary  to  characterize  these  coefficients  and  test  the  model. 

6.  CONCLUSIONS 

In  summary,  we  have  measured  the  excited  state  properties  of  donor-acceptor  push-pull  charge-transfer  chromophores 
and  modeled  the  nanosecond  NLT  as  an  effective  three-photon  absorption  process  combined  with  stimulated  scattering 
at  high  intensities.  We  have  also  considered  the  effects  due  to  nonlinear  refraction  by  employing  a  three-dimensional 
nonlinear  beam  propagation  model.  Several  model  parameters,  including  intrinsic  TPA  cross  sections,  have  been  meas¬ 
ured  independently,  and  the  analytical  model  shows  relatively  good  agreement  with  NLT  data  from  solutions  for  which 
the  molecular  TPA  and  ESA  cross  sections  have  been  accurately  measured.  Discrepancies  appear  when  the  excited 
singlet-singlet  cross  sections  are  not  well  determined  (i.e.,  the  excited  state  absorbance  measurements  are  noisy),  and 
smaller  discrepancies  are  also  observed  in  almost  every  sample,  particularly  in  the  low  energy  region.  We  have  argued 
that  it  is  unlikely  that  any  of  these  discrepancies  originate  from  stimulated  scattering  effects,  and  that  it  is  plausible  that 
their  origin  may  rooted  in  nonlinear  refraction.  We  conclude  that  the  dominant  contribution  to  the  nonlinear  transmis¬ 
sion  in  these  chromophores  in  the  nanosecond  regime  is  ESA  from  singlet  states,  with  only  minor  contributions  from 
triplet  states  due  to  the  low  triplet  quantum  yield.  TPA  does  not  contribute  significantly  to  the  transmittance  loss,  but  is 
key  to  pumping  the  excited  states.  Nonlinear  refraction  must  be  studied  in  more  detail  to  determine  its  impact  on  the 
NLT  at  both  low  and  high  intensities.  Further  modeling  work  is  required  to  properly  incorporate  nonlinear  refraction 
from  population  redistribution  and  thermal  effects,  and  to  measure  these  coefficients  in  sample  solutions. 
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